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ABSTRACT 

Our  aim  in  this  work  is  to  show  that,  in  a  "permanent  regime",  the 
behaviour  of  a  viscous  incompressible  fluid  can  be,  in  principle,  determined 
by  the  study  of  a  finite  number  of  modes.  It  is  proved  that  the  behaviour 
for  t  -*■  “  of  the  solution  to  the  Navier-Stokes  equations  is  completely 
determined  by  its  projection  on  appropriate  finite  dimensional  subspaces, 
corresponding  to  eigenspaces  of  the  linear  operator,  or  more  general 
subspaces,  including  finite  element  subspaces.  Some  indications  on  the 
dimension  of  such  subspaces  are  given. 
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SIGNIFICANCE  AND  EXPLANATION 


If  a  viscous  incompressible  fluid  is  driven  by  time  independent  forces  of 
sufficient  intensity  then,  after  a  transient  period,  the  "permanent"  regime 
seems  to  be  totally  chaotic  and  unstructured.  The  present  work  is  part  of  a 
set  of  articles  which,  however,  tend  to  show  that  there  may  be  some  structure 
in  such  flows  (at  least  in  the  case  where  the  space  dimension  is  2),  in 
particular  that  they  are  determined  by  a  finite  number  of  parameters. 
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A^^MPTOTIC  NUMKKICAL  ANALYSIS  FOR  THE  NAVIER-STOKKS  EQUATIONS  (I) 


C.  Foias  and  R.  Temam 


Introduction 

Up  to  now  the  numerical  analysis  or  the  Navier-Stokes  equations  has  been  limited  to 
tne  study  ot  the  approximation  of  time  dependent  solutions  on  a  finite  interval  of  time  ot 
to  the  approximation  of  stationary  solutions  (cf.  among  many  other  references  [1,4, 8, 9]). 

In  the  presence  of  a  turbulent  flow  driven  by  a  steady  excitation  a  different  type  of 
problem  arises  naturally:  the  study  of  the  long  time  behavior  of  the  solutions. 

The  present  is  an  essay,  the  purpose  of  which  is  to  show  that  for  the  2-D  Navier- 
Stokes  equations  and,  under  some  circumstances  for  the  3-D  Navier-Stokes  equations,  there 
is  a  theoretical  basis  for  determining  the  qualitative  long  time  behavior  of  a  fluid  by  the 
study  of  a  finite  number  of  adequate  modes,  A  typical  result  is  the  following  one:  Let 
W  be  a  finite  dimensional  subspace  of  the  natural  function  space  V.  If  W  satisfies  a 
certain  condition,  then  the  behavior  for  t  +  “  of  a  solution  u  of  the  Navier-Stokes 
equation  is  completely  determined  by  the  behavior  for  t  +  ”  of  its  projection  on  w. 

Several  results  of  this  type  are  derived  in  this  article.  While  this  kind  of  problem 
was  already  discussed  in  [21,  our  present  interest  was  aroused  and  inspired  by  the 
questions,  conjectures  and  ideas  due  to  0,  P,  Manley  and  y,  M.  Treve  [2,12,11),  with  whom 
we  acknowledge  fruitful  discussions  and  correspondence.  In  this  paper  we  did  not  try  to 
produce  the  best  constants,  and  we  did  not  try  to  present  the  main  inequalities  in  a 
nondimensional  form.  These  improvements  of  the  work,  and  other  developments  will  appear  in 
a  subsequent  work  [13], 

The  plan  is  as  follows: 

1,  Notations  and  recapitulation  of  results, 

2,  Approximation  in  the  subspaces 

3,  Approximation  in  a  general  subspace. 

4,  Time  periodic  solutions. 

8,  Remark  on  Galerkin  approximation. 


Sponsored  by  the  United  States  Army  under  Contract  Ho.  DAAG29-80-C-P041 . 


1 .  Notations  and  recapitulation  of  results 

Let  it  be  a  bounded  domain  of  R  ,  t  =  2  or  3,  and  let  T  be  its  boundary.  We 
assume  that  f  is  a  manifold  of  dimension  4-1,  of  class  with  a  finite  number  of 

connected  components,  and  that  0  is  locally  located  on  one  side  of  f.  We  shall  firstly 
consider  the  initial  value  problem  for  the  tJavier-Stokes  equations: 


(1.11 

(1.i) 

(1.3) 

(1.4) 


-  viu  +  (u*7)u  +  Vp  »  f  in  £1 
7u  =  0  in  (1 
u  =  0  in  r 


“'t-O  ■  “0  ' 

where  v  >  0  is  the  (cinematic  viscosity,  u  =  (u^,U2)  or  (u^,U2,Uj)  is  the  velocity,  a 
vector-valued  function  of  x  G  0  and  t  >  0,  and  p  is  the  pressure,  p  =  p(x,t),  and 
f  represents  the  external  laody  force  per  unit  of  mass. 

All  what  follows  apply  to  the  case  where  (1.3)  is  replaced  by  a  nonhomogeneous 
boundary  condition,  which  corresponds  to  more  realistic  physical  situations  (Couette-Taylor 

flow,  (91).  We  will  refrain  fror  treating  this  case  to  avoid  purely  technical  difficulties. 
2 

We  denote  by  L  (fl)  the  space  of  square  integrable  real  functions  on  S  and  by 

1  2 
Hy(fl)  the  Sobolev  space  made  of  the  funetions  which  are  in  L  (0)  together  with  their 

2  2  t 

first  derivative  and  which  vanish  on  T,  we  set,  for  u  in  L  (£1)  or  L  (£))  (resp. 
u  in  H^(S))  or  HQ(il)^) 

|u|^  =  /  lu(x)|^dx,  resp.  lul^  =  /  |Vu(x)|^dx  , 

n  £2 


2  4 

The  space  L  (£1)  admit  classically  an  orthogonal  decomposition  of  the  form  H  ®  G,  where 

G  =  {v  »  Vq,  q  G  h'(£1)} 

and  its  orthogonal  H  satisfies 

H  =  {v  G  L^(£l)*,  div  v  =  0,  v*n|j.  =  0}  , 

n  the  unit  outward  normal  of  F.  Let  also 

1  4 

V  =  (v  6  H|j(£l)  ,  div  V  »  O)  . 

It  IS  clear  that  H  and  V  are  Hilbert  spaces  for  the  norms  1*1  and  •*•,  and  their 
corresponding  scalar  products.  While  |v|^  is  equal  to  the  Itinetic  energy  of  the  fluid 
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with  velocity  V  (the  density  P  =  1),  we  recall  that  for  v  e  V,  Ivl  too,  reduces  to  a 
physical  quantity: 

«v«^  =  /  [curl  v(x)l^dx  , 

n 


2  i 

Let  It  denote  the  orthogonal  projection  of  L  (SI)  onto  H  and  define  the 
operators  A  and  B  by 

(1.5)  Au  =  -itAu  for  u  G  0(A)  .  V  n  H^(n)*  , 

(1.6)  B(u,v)  c  tr((u,7)v),  for  u,v  e  D(A)  . 

Then  A  is  a  self-adjoint  operator  in  H  with  an  orthonormal  basis  jw  }  of 

m>1 

eigenvectors,  such  that 

(1.7)  Aw„  =  X^w_^,  m  >  1,  0  <  <  Xj  <  ...  , 


X  as  m  ♦  +“  and 

m 

(1.8)  V  =  D(A  ^2  ),  lul  «  |a  <^201  for  u  6  V  . 

The  operator  A  is  an  isomorphism  from  D(A)  onto  H  and  from  V  onto  V’  (the  dual 
of  V  which  one  can  identify  to  a  superspace  of  H).  Concerning  B  we  recall  the 
following  fact:  B  is  a  compact  mapping  from  D(A)  x  V  or  V  x  D(A)  into  H  and  from 
V  X  V  into  V.  Furthermore  we  have  the  estimates: 


(1.9) 


(1.10) 


lB(u,v)»^,  < 


/2  iur/^«ul’/2|vp/^v«’''2 


B(u, v) I  < 


|ul’''^lAut^^^»vl 

|u|'^2|ui’/2,^,V2.  ,1/2 


if  I 
if  f 


2 

3 


if  1=2  and  if  f  =  2  3: 

f  Iau]  KvI 

(1.11)  lB(u,v)j  <  c,  < 

[  IuHIavI  . 

The  constants  Cq  and  c,  depend  only  on  (2  but  are  not  easy  to  determine  since  they 
invo  ve  the  norm  of  the  operator  A”',  and  Sobolev  constants.  All  the  above  results  can 
be  found  for  instance  in  [6J,  [91. 

We  can  now  recall  the  functional  formulation  of  (1.1)-(1.4):  This  is  the  functional 
differential  equation 
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(1.12)  ^  +  VAu  +  B(u,u)  =  £,  “It^O  “  “O  ' 
where  we  shall  assume  for  simplicity  Uq  €  V  and 

£  is  continuous  and  bounded  from  [0,“)  into  H, 

(1.13)  •  2 

£'  =  —  belongs  to  L  (0,“;V')  . 
dt  loc 


It  is  well-)tnown  that,  if  4  =  2,  there  exists  a  unique  function  u  such  that 

(1.14)  ue  C(  IO,“);H)  n  L^^^(  tO,“);V) 

satisfying  the  equation  (1,12)  in  V*  (cf.  for  instance  [9]).  Moreover  u  is  actually  a 
continuous  function  from  (0,*)  into  D(A),  which  is  bounded  in  D(A)  as  t  +*•  We 
set 

(1.15)  Bu(t)i  <  C2  “  ^  t  >  0  , 

(1.16)  |Au(t)|  <  C3  =  C2(|u^|,lfl,1/v,C|j,X^,1/a)  for  t>a>0 
where 

(1.17)  [fl  =  sup  |f(t)|  . 

0<t<" 


The  estimate  (1.15)  is  given  in  (21;  the  estimate  (1.16)  is  more  recent  and  given  in  151 
(cf .  also  [101 ). 

if  4  =  3,  there  exists  a  wealt  solution  of  (1.12)  bounded  in  H.  Such  a  solution  may 
or  may  not  be  bounded  in  V.  We  will  only  consider  such  a  solution  if  it  is  bounded  in  V 

(1.18)  lu(t)l  <  R  for  t  >  0  . 

In  this  case  it  follows  also  from  [31,  151,  110),  that  u(t>  belongs  to  D(A)  and  is 
bounded  in  D(A)  as  t  ♦  <“: 

(1.19)  |Au(t)|  <  c^(R,lf),1/v,c^,X^,1/o)  for  t  >  o>  0  . 

In  the  sequel  we  will  let  £  *  2  or  3,  and  consider  solutions  of  (1.12)  which  are 
uniformly  bounded  in  V  on  (0,*);  the  existence  of  R  <  •  is  an  assumption  if  i  =*  3, 
and  is  automatic  if  i  *  2  (R  »  C2  cf.  (1.15)  )•  See  the  comment  in  Remark  2.1. 
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2.  Approximation  In  the  aubspaces 
2.1.  An  inequality 

For  m  e  N,  we  denote  by  the  space  spanned  by  the  eigenfunctions  w^,..,,W|^, 

of  A,  and  by  Pjjj  the  projector  in  V,  H  or  V  onto  Finally  =  I  -  It  is 

easy  to  see  that 

(2.1)  I ^  jv>  I .  e 


(2.2)  Itl  <  V*  6  tt^V  . 


Let  u(*) 

and  v(*)  be  solutions  of 

the 

equation 

(2.3) 

u'  +  VAu  +  B(u,u)  « 

f. 

t  >  0, 

u(0)  »  u^ 

(2.4) 

V'  +  UAv  +  B(v,v)  “ 

9/ 

t  >  0, 

v(0)  •  V. 

where  6  V  and  g  satisfies  the  same  assumption  (1.13)  as  f.  We  set 

w  »  U  -  V,  p_  »  P_w,  q_  »  Q_w,  e  »  f  -  g,  e_  =  O  e  . 

m  '  *in  *  ’ '  m  *in 

Then 


and  consequently* '  \ 


(2.5) 


-  +  (e„,q^)  . 

The  right-hand  aide  of  (2.5)  can  be  bounded  because  of  (1.11)  by*^*! 


*’*we  recall  that  (B(  >^,4i)  ,9)  •  -(B(  <>,  0) ,  i|i)  for  in  V. 

( 2 ) 

Ci»c|,c|[, . . . ,  denote  various  positive  constants.  The  c^'s  are  the  same  all  the  time, 
may  represent  different  quantities  at  different  places. 


i 
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<  (for  E  >  0  arbitrary) 

,-1  2 


^  r  'V'  ^ ^  ^ f  'V'  * ^  ^  ^  '^iCf 


<  (by  (1.16),  (1.19)) 


-1 


-  2  2 
2c.  c. 


<  EvJq  1^  +  — j—  |e  1^  +  — ^  Ipml^  +  c  c  X  '{^Iq  1^ 
^  2ev  '  m'  e\)  1  4  m+1  m 

for  t  >  a  , 

We  have  thus  proved  the  following 
Lemma  2 . 1 


If  m  is  large  enough,  so  that 


(1) 


(2.6) 


2  2 

‘=1^4 
X  >  '  * 

m+1  *  .2  ' 


then,  for  t  >  a  >  0, 


(2.7) 


,'1 


^22 
4c. c , 


l<lml"  -  ‘"'V'  *  if  l«ml'  -  IPml^ 


c  c  X'’/^ 

y  =  V  -  c,o  X-V2  ,  0,  c=i--i4J2il- 

1  4  m+1  2  2V 

2.2.  The  main  result 


We  introduce  the  following  weak  mode  of  convergence:  we  will  say  that  ‘^(t) 
converges  essentially  to  0  as  t  ♦  “  (and  we  write  i^(t)  l  +  0)  if: 

(2.8)  there  exists  x  >  0,  such  that  for  every  5  >  0,  there  exists 

tj  satisfying 

measCi  6  (t,t  +  l),|y>(t)|  >  6}  <  <5,  for  every  t  >  tj  . 


(1  ) 


36c. 


We  need 
2 


^  ,  * 
m+1 


(sup  |Au(t)t^) 
t>o 


X  ^  ‘  ®“P  |hu(t)|^. 


t*** 


and  actually  it  is  sufficient  to  have 
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It  is  easy  to  see  that  if  ^(t)  ♦  0  for  t  +  then  ^(t)  o,  for  t  +  *.  Also 

if  V?  e  l'(0,®)  Chen  ^(t)  —■—*■•>  0  as  t  ♦  *. 

Theorem  2.1 

We  assume  that  1=2  or  that  i  «  3  and  that  u  and  v  are  solutions  of  (2,3), 
(2.4)  uniformly  bounded  in  V.  We  assume  also  that  (2.6)  is  satisfied.  Then; 
i£  "  V(t)|  +  0,  |(I  -  P^)(f(t)  -  g(t))|  +  0,  for  t  +  »,  then 

(2.9)  |(I  -  p^)(u(t)  -  v(t))(  ♦  0  for  t  ♦  * 

(2.10)  |u(t)  -  v(t)l  ♦  0  for  t  . 

ii)  21  IPm'i'lt)  ■  v(t))l^  0,  j(I  -  P|^)(f(t)  -  g(t))l^  0,  for  t  +  then 

(2.9)  still  holds  but  instead  of  (2.10), 

(2.11)  |u(t)  -  v(t)  0  for  t  ♦  -  . 

Proof 

We  infer  from  (2.2)  and  (2.7) 


d_ 

dt 


+  VX’/f|q  ! 
m+1  '^m 


2  2 
4c  c 

CV 


iPml  • 


Whence  for  t  >  t^  >  o; 


(2.12) 


h„(t)|  <  lq„(to)|  e 


.,2  2 

^^1^=4 


ev  '‘^ro 


|p„(T)|2]e 


dT 


There  exists  k,k'  >  0  such  that  for  every  6  >  o,  there  exists  t^  which  satisfies 
for  t  >  tj  : 

measCi  e  (t,t  +  1),  |p  (u(T)  -  v(T))l^  >  S}  <  k6 

mea8{T  e  (t,t  +  1),  |p^(f(T)  -  g(T))|^  >  6}  <  k'5  . 

For  a  fixed  integer  M  we  ta)te  t  >  tg  +  H,  tg  >  max(tj,a).  Then  (2.12)  implies 
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X"’  4c^o^  t 

t-M 


<  c(u,v)e  ”’■*■’  “  +  «(^ 

X"’  4c?c?  t-M  -Vx'''^?(t-T) 


+  c(f,g)  +  **  c(u,v)]  •  (/  e 

*^0 


m+1 


dTj 


, -1 


*  2  2 
4c. c  , 


+  m6(k'  — —  c(f,q)  +  <  — c(u,v)] 


where 


c{u,v)  =  supjuCt)  -  v(t)l,  c(£,g)  =  8up|f(t)  -  g(t)| 
t>0  t>0 


Therefore,  as  t  ♦ 


2  ^^4  1 

Um  sup|q  (t)p  <  4(—  +  -^)( - -^72 - ) 

t+«  v’X  ' . 

Bl+1 


,  -1 


X-'  4c2c2 

[^c(f,g)  +  — —  c(u,v)](^ - ^-7^) 


v.xV2 

m+1 


v-1 


-22 

4c, 


+  m5[k'  c(f,g)  +  <  — c(u,v)]  . 


EV  V 

We  let  6*0  and  then  M  ♦  and  we  obtain  (2.9). 


Remark  2.1 


This  result  which  is  contained  in  a  slightly  weaker 
space  is  4=2,  is  reproduced  here  for  the  convenience 


form  in  (21,  when  the  dimension  of 
of  the  reader. 
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3.  Approximation  in  a  general  subspace 

In  the  applications  the  utilization  of  the  basis  [w  }  formed  by  the  eigenvectors 

of  A  13  not  practical  since  these  functions  are  not  easy  to  determine.  Therefore  we 
shall  now  show  how  the  preceding  theorem  can  be  extended  when  is  replaced  by  a  general 

finite  dimensional  subspace  W  of  V. 

3.1.  Assumptions 

Let  W  be  a  finite  dimensional  subspaee  of  v  and  let  P(W)  be  the  pro]ector  in 
H  onto  W,Q(W)  =  I  -  P(W).  Since  P(V))  is  not  a  projector  in  V,  it  may  happen  that 
((t'/'P))  ^0  if  ^  S  W  and  iji  6  V,  P(W)<ji  »  0.  However,  one  can  show  (cf.  Lemma  3.2 
below)  that  there  exists  P(W),  0  <  p(W)  <  1  such  that 

(3.1)  I  ( ('<>,4') )  I  <  p(W)«v>»#4'«,  vv>  e  w,  vji  e  v,  p(w)4i  =  o. 

We  associate  also  to  W  the  two  numbers  A(W),p(W) 

X(W)  =»  inf{Sv'i^,4’  s  V,  P(W)v!  =  0,  |v;|  =  1} 

U(W)  =  SupO'CI^,  S  W,  ||(||  =  i)  , 

so  that 

(3.2)  |>p1  <  VV>  e  V,  P(W)^  =  0 

(3.3)  «  P(W)’^^|>|il,  V<|i  6  W  . 

When  it  is  not  necessary  to  mention  the  dependence  on  w,  we  will  write  simply 
P,Q,P/X,U,  instead  of  P(W),...  . 

3.2.  An  inequality 

We  consider  as  in  Section  2,  the  solutions  u<*),v(')  of  (2.3),  (2.4)  and  we  set 
w  =  u  -  y,  p  =  pw,  q  =  Qw,  e  =  f  -  g  . 

We  ha  ve 

q'  +  vQAw  +  QB(v,w)  +  2B(w,u)  =  Qe 
and  talking  the  scalar  product  in  K  with  q,  we  get 

^  |q|^  +  Vllqll^  =  -(B(v,p),q)  -  (B(p,u),q), 

-(B(q,u),q)  +  (Qe,q)  -  V((p,q))  . 

Using  (1.11)  and  (3.1)-(3.3)  we  find  that  the  right-hand  side  of  this  inequality  is  less 
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than 


XJge  I  Hql  +  |p1  Iql  +  (  1au|  +  |av|  )  |  p|  Iql  +  c  Jau|  jqj  »q  I 


«  (by  (1.16),  (1.19)) 


,  e\<  .  ,2  1  ,  .2  ev  .  ,2  3V  2  ,  .2 

1“ +  — 'q'  +  —  p  kIpI 


ic 

>  -p  Iq*^  +  (|au|  +  |av|)^|p|^  + 


,,-1  2  2 
.2  ^^1  ,2  .  3  ,‘^l‘^4  .  ..2..,,  ,2 


<  £vlql  +  |i)e|  *  -  (-^  +  yp  p)|p|  +  c^c^X"  Iql  tor  t  >  a  >  0  . 


tfien  we  set 


A(W)  > 


1/2  1  C  C  X(W)"’'^^ 

y  =  y  -  c,c^X(w)-  /  >0,  e=- - - - 


and  we  have  established: 


Lemma  3. 1 


If  (3.1),  (3.4)  hold,  then  for  t  >  o  >  0, 


^  |q|2  +  ylq|2  <  ^  jgel^  +  I  +  vp2u)|p|2  , 


y  ,  e  as  in  (3.5). 


3. 3.  statement  of  the  result 


As  in  Section  2.2,  we  have 


Theorem  3.1 


We  assume  that  t  =  2  or  that  t  =3  and  that  u  and  v  are  solution  of  (2.3), 

(2.4)  uniformly  bounded  in  V.  Let  W  be  a  finite  dimensional  subspace  of  V  such  that 

(3.4)  is  satisfied.  Then 

i)  |?(u(t)  -  v(t))l  -►  0,  1(1  -  P)(f(t)  -  g(t))|  +  0  for  t  ♦  then 

(3.7)  |(I  -  P)(u(t)  -  v(t))l  *  0  for  t  +  ®  , 

(3.U)  lu(t)  -  v(t)|  ♦  0  for  t  +  "  . 

11)  JJ_  |p(u(t)-v(t))|^  0,  1(1  -  P)(f(t)  -  g(t))l^  ■— *°*-*  0  for  t  ♦  ",  then 


(3,7)  holds  .^nd,  instead  of  (3.8) 
(3.9) 


u(t)  -  v(t)l  - 0,  .for  t 


The  proof,  stsrtinij  tram  (3.b),  is  essentially  the  same  as  that  of  Theorem  2.1. 

In  the  rest  of  this  section  we  give  examples  and  show  that  (3.1)  is  always  satisfied, 
while  assumption  (3.4)  is  satisfied  if  w  is  “sufficiently  large". 

3.4.  Assumption  (3.1). 

Lemma  3.2 

Under  the  assumptions  of  Section  3.1,  there  exists  p  =  P(w),  0  4  p  <  1  such  that 
( 3.1 )  holds. ' 

Proof 

If  (3.1)  was  not  true,  we  could  find  two  sequences  (.f.  }  #  I't'- }  >  6  W,  1(1.  6  V, 

3  j>i  ^  j>1  J  ^ 

=  0,  such  that 

Ivi.lllili.i  >  !((■<’.,</.))  I  >  (1  -  4-)I^.  I  . 

3  3  3  3  3  3  3 


Setting  we  find 

3  3 


(3.10) 


1  >  |<(4>',i|-'))i  >  (1  -  j)  . 


We  can  extract  a  subsequence  (still  denoted  j)  such  that  4>j  converges  to  some  limit 
ll^H  =  1 ,  .s  6  w  (w  has  finite  dimension),  and  i|)j  converge  wea)tly  in  V  to  ♦,  iji  e  V, 
li)/l  <  1  ,  pi(i  =  0.  At  the  limit,  (3.10)  gives 

|((^,i)'))|  =  1,  ■=  1,  1*1  4  1  , 

so  that  =  1,  =  ici<:  *  0,  by  contradiction  with  P*  =  0. 

3.5.  Assumption  (3.4)  -  Lxample 

We  consider  the  following  situation  which  is  classical  in  the  numerical  analysis  of 
partial  differential  equations  and  in  particular  of  Navter-Sto)ces  equations  (cf.  [91,  Chap. 
I,  ^4): 


(  1  ) 

any  role  in  t)ie  proof  above.  The  inequality  (3.1)  with  p(W)  =  1  is  trivial. 


Actually  the  tact  that  P  <  1  which  is  important  in  other  developments,  did  not  play 
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We  are  given  a  family  (w  }  i  of  finite  dimensional  subspaces  of  V. 

^  heH 


The  set  of 


indices  H  is  arbitrar’^  but  is  equiped  with  a  concept  of  limit* which  we  denote  for 

simplicity  lim.  For  example,  in  the  Galerkin  method,  H  =  l/Mi  h  =  V®,  and  we  pass  to 
htO 


the  limit  m  +  h  ♦  0,  For  finite  element  methods  (cf.  191),  H  is  a  family  of  regular 
triangulations  of  the  domain  fl),  and  we  let  the  diameter  of  the  largest  triangle  go  to  0. 

The  main  assumption  on  the  spaces  Wj^  is  the  following  one 
(3.11)  WVJ  6  V,  Int  as  h*0. 


^h 


In  the  case  of  an  increasing  sequence  of  subspaces  of  V  (Galer)cin  method), 

assumption  (3.11)  means  simply  that 

(3.12)  LJ  W  is  dense  in  V  . 

in€M  " 


Then  we  see  that  assumption  (3.4)  is  satisfied  for  h  ‘sufficiently  small". 
Lemma  3.3 

Under  assumption  (3.11), 

(3.13)  lim  X(W^)  -  +«  , 

h*0 

and  (3.4)  is  satisfied  for  h  sufficiently  small. 

Proof 

The  proof  consists  in  showing  the  more  precise  following  statement 

(3.14)  For  every  integer  m,  there  exists  h„  and,  for  h  <  h  ,  A(w.  )  >1  -  1. 

m  n  m 

For  given  m  and  6  >  o,  the  assumption  (3.11)  written  with  ip  =  w^,  j  = 

shows  that  there  exists  h.  such  that 

m 

Inf  Iw.  -  <  6,  tor  j  *  and  for  every  h  <  h  • 

1  '  m 

h 


*’ ’a  filter  F 


with  a  denumerable  basis; 


lim  means  roughly  spealting  lim. 
h»0  F 
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The  notations  are  the  same  as  in  Section  3*  If  the  assumption  (3.4)  is  satisfied,  we 
infer  from  Lemma  3.1  that  there  exist  two  positive  constants,  >  0,  n  >  0,  independent 

of  u  and  v  such  that 

-n(t-t  ) 

(4.1)  ({I  -  P)(u(t)  -  v(t))|^  <  ^5^® 

+  J  e"''“^"^*l|p(u(T)  -  v(T)j^  +  |(I  -  P)(f(T)  -  g(T))|^]dT)  , 

"o 

for  t>t|^>a>0. 

We  can  now  prove  the  following: 

Theorem  4. 1 

We  assume  that  ^  =  2  or  that  i  =  3  and  that  u  is  a  solution  of  (2.3)  uniformly 
bounded  in  V.  Let  W  be  a  finite  dimensional  subspace  of  V  such  that  (3.4)  is 
satisfied.  Assume  moreover  that  there  exist  periodic  functions  f»(*)>  with  value 

in  H  and  W  and  period  T  >  0,  such  that 

(4.2)  lim  lf<t)  -  fjt)|  =  lim  |Pu(t)  -  p„(t)|  =  0  . 

t-H» 

Then  there  exists  a  periodic  solution  u^  with  period  T  of  the  equation 

(4.3)  “1  +  +  B(u  .u  )  =  f  , 

09  09  09  09  09 

such  that 

(4.4)  lim  iu(tj  -  u^(t)l  =  0  , 

Proof 


Apply  (4.1)  to  u  and  v,v(t)  »  u(t  +  jT).  Then  for  JtT  >  tQ, 


we  obtain  at  time 


Id  -  P)(u(t  +  it)  -  u(t  +  (J  +  »)T))|  < 


<  c^{e 


+  i  e 


•  l|P(u(T)  -  u{T  +  jT))|^  +  1(1  -  P)(£(T)  -  £(T  +  jT))|^]dT} 
For  e  >  0  given,  let  >  tg  be  such  that  for  t  >  t^ 

IPvi(t)  -  P„{t)l  <  e,  |f(t)  -  f„(t)l  <  E  , 

Then  from  (4.5)  we  obtain  that 


-n(t+fT-t  ) 

|(I  -  P)<u(t  +  fT)  -  u(t  +  (j  +  f)T))|  <  CjB 


'"C  ^4-8.T 

/  I  -n(t+iT-T)  2  f  -n(t+*T-T). 

+  CjcCu.f)  J  e  dT  +  8c^e  J  e  c 


2  *nfi.T 

sup  1(1  -  P)(u(t  +  tT)  -  u(t+(4  +  j)T))l  <ce  +  — , 

tg 

for  every  j  *  0  and  i  >  — ,  Therefore  {u(t  +  W) )  is  a  Cauchy  sequence  in  the  space 
of  continuous  founded  functions  from  (tg,“)  into  (I  -  P)W,  Thus  there  exists  a 
continuous  bounded  function  u  from  (t„,<*)  into  H,  such  that 

«0  Q 

(4.6)  u(t  +  It)  ♦  U|^(t)  as  I  ♦  in  H,  uniformly  in  t  on  [tg,“)  . 

Since  lAu(t  +  iT)  1  <  c,j  for  all  t  >  a,  f  >  1 ,  we  see  that  |Au^(t)l  <  c^  for  all 


t  >  a,  and 


u(t  +  iT)  ♦  u  (t)  in  V,  V,  t  >  a  . 


It  is  then  easy  to  see  that  u^  is  a  solution,  bounded  in  v,  of  (4.3),  and  u^  is 

periodic  of  period  T,  just  because  of  (4.6).  The  convergence 

(4.8)  lim  lu(t)  -  u^(t)i  »  0  , 

t*» 


follows  immediately  from  (4.6)  as  well  as  the  T  periodicity  of  u^(*).  i  r^e  convergence 
(4.4)  follows  from  (4.8)  and  the  fact  that  |Au(t)  -  Au^(  remains  bounded  (by  2c^)  as 
t  ■»  «.  ■ 

We  then  deduce  the  following  result  for  stationary  solutions 
Theorem  4.2 

The  assumptions  are  similar  to  that  of  Theorem  4.1.  We  assume  that  there  exist 
f^  6  H  and  p^  6  W  such  that 

(4.9)  lf(t)  -  f^l  0,  |pu(t)  -  p^l  ♦  0,  for  t  +  ••  . 

It  follows  that  there  exists  e  D(A)  such  that 

(4.10)  lu(t)  -  u^l  0  for  t  ♦  “  , 
where  u^  is  a  (stationary)  solution  of  the  Navier-Stolces  equation 

(4.11)  VAu^  +  B(u^,u_)  =  f_  . 

Ok'  M  a» 

Proof 

We  apply  Theorem  4.1  with  T  >  0  fixed,  arbitrary,  and  we  obtain  (4.3),  (4.4).  Now 
T  >  0  can  be  chosen  arbitrarily  small  and  since  u^  must  be  independent  of  T  in  (4.4), 
we  conclude  that  u_^  has  period  0,  i.e.  u^  is  independent  of  t.  • 


-16- 


Remark  on  Galerkin  approxi mation 


We  now  assume  that  the  dimension  i  2. 

For  simplicity  we  restrict  ourself  to  a  Galerkin  approximation  of  Navier-Stokes 
equations  based  on  the  spaces  V^,  i.e.  the  family  of  eigenfunctions  of  the  Stokes 

problem.  We  will  show  that  if  m  is  sufficiently  large,  the  behavior  as  t  *  “  of  the 
Galerkin  approximation  Uj^^,  is  completely  determined  by  the  behavior  as  t  +  *  of  a 
certain  number  m,  of  its  modes,  i.e.  of  Pjjj  u^^,  m,  <  m.  This  number  m«  does  not  depend 
on  m. 

5.1.  Galerkin  approximation 

For  fixed  m,  the  Galerkin  approximation  Uj^^  of  the  solution  u  of  (2.3)  is  defined 

bys 

f  “m  ^  ''*"m  V<“m'“m>  “  V<  ^  °  - 

1  "m<0>  ■  Vo  • 

It  is  classical  to  derive  a  priori  estimates  Independent  on  m  on  u^^:  for  example, 
for  every  t  >  s  >  0: 


(S.2) 


,  -  A,  (t  -  s) 

lu^(o)rdo  <  l%<s)r  - y- -  I£1 


and 


(5.3) 


-vX,t 
e  ) 


(f  )■ 


The  following  a  priori  estimate  is  verified  by  Uj^: 

Lemma  5.1 

•u^(t)l  is  bounded  independently  of  m  and  t  for  t  >  a  >  0,  m  >  0. 

Proof 


Taking  the  scalar  product  of  (5,1)  with  we  obtain 


T  ^  +  v|au 

2  dt  m  '  m' 


'  ♦tt.Au^)  -  (B(u„,u„ ) ,Au„  . 
m  m'  m  '  m 


Because  of  (1.10),  (1.17),  the  right-hand  side  of  this  expression  is  majorized  by 
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|Q_  6||q  I  +  c'|v  I’^^lq  I 

*  m'  m  ’  nt. 


«*  •  ‘  «»* 


+  C'|v  I’/^lq  I’^^lq  I^'^^+c'Ip  « Iq^  I 


m.  m.  m 


If 

(5.6) 
we  set 


(5.7) 


c-r’/'* 

\.<  -  2(V  -  C“X'’{^)  >0,  e  -  1 - >  0  , 

W^+l  4i  <2v 


and  we  bound  the  last  quantity  above  (5.6)  by 


,-1 


ev  „_  „2  .  1 


(c“)^xy^ 

2  .  .-  ,2  .  ----*-  Ip  .  cn-VSq  ,2 

2CV  '  *^v<.  m  4.1  ^WI 


and  we  find 


(5.B) 


|2^v.|q  |2«-i|Q  e|2.^  - !lL|p 

dt  '  ev  '''m.  '  EV 


As  for  Theorem  2.1  we  obtain 
Theorem  5.1 

We  assume  that  4  «  2  and  that  m  >  m,,  ro,  sufficiently  large  so  that  (5.6) 
verified.  Then: 
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i)  If 


,(t))|  *  0.  |(I  -  )(f(t)  -  g(t))|  ♦  0 


for  t  +  »»;  then 

(5.9) 

(5.10) 

ii)  If 


|(I  -  ~  ^  t  +  <•  , 

|u^(t)  -  V||,(t)l  +0  for  t  ♦  ®  . 

-  v^(t)) 0,  |(I  -  P^^)(f(t)  -  g(t))|2  0 


for  t  ♦  <”,  then  (5.9)  holds  and  instead  of  (5.10), 

(5.11)  |u_^(t)  -  v^(t)|2  0,  for  t*-. 

Remar)c  5.1 

This  theorem  will  ta)te  its  full  Interest  if  we  can  relate  the  behavior  for  t  +  “  of 
the  Galer)tin  approximation  of  u,  to  the  behavior  for  t  “  of  u  itself.  This 

question  will  be  considered  in  a  subsequent  work. 
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